Spin Susceptibility in Underdoped YBa2Cu306n 
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We report a comprehensive polarized and unpolarized neutron scattering study of the evolution of 
the dynamical spin susceptibility with temperature and doping in three underdoped single crystals 
of the YBa2Cu30e+ x high temperature superconductor: YBa2Cu30e.5 (T c = 52 K), YBa2Cu306.7 
(T c = 67 K), and YBa2Cu30e.85 (T c = 87 K). The spin susceptibility is determined in absolute units 
at excitation energies between 1 and 140 meV and temperatures between 1.5 and 300 K. Polarization 
analysis is used extensively at low energies. Transitional matrix elements, including those between 
spin states, in a bilayer system such as YBa2Cu306+ x can be generally classified into even and odd, 
according to the sign change under a symmetry operation that exchanges the layers, and both even 
and odd excitations are detected in YBa2Cu3C>6.5 and YBa2Cu306.7. While the even spin excitations 
show a true gap which depends on doping, the odd spectrum is characterized by a weakly doping- 
dependent pseudogap. Both even and odd components are substantially enhanced upon lowering 
the temperature from 300 K. The even excitations evolve smoothly through the superconducting 
transition temperature T c , but the odd excitations develop a true gap below T c . At the same 
time, the odd spin susceptibility is sharply enhanced below T c around an energy that increases with 
doping. This anomaly in the magnetic spectrum is closely related to the magnetic resonance peak 
that appears at 40 meV in the superconducting state of the optimally doped compound (T c = 93 K). 
From these data we extract the energy and the energy-integrated spectral weight of the resonance 
peak in absolute units as a function of doping level. Theoretical implications of these measurements 
are discussed, and a critique of recent attempts to relate the spin excitations to the thermodynamics 
of high temperature superconductors is given. 



PACS numbers: 74.25.Ha, 74.60.W, 25.40. Fq, 74.72.Bk 



I. INTRODUCTION 

The importance of electronic correlations in the cop- 
per oxide high temperature superconductors is now 
generally recognized. Some of the strongest evidence 
has been provided by neutron scattering measurements 
of the spin excitations in several families of cuprates: 
La 2 _ x Sr x Cu0 4 0, YBa 2 Cu 3 6+x (2j-[2TJ , and recently 
also Bi 2 Sr 2 CaCu208+5 J22|. While in an uncorrelated 
metal the spin excitation spectrum takes the form of a 
broad continuum extending up to energies comparable to 
the Fermi energy, the spin excitations in the metallic cop- 
per oxides in many ways resemble the antiferromagnetic 
magnons in their insulating parent compounds. This is 
illustrated by the mere fact that spin excitations cannot 
be observed at all in the metallic materials by neutron 
scattering: the spectral weight of continuum excitations 
in an uncorrelated metal is more than an order of mag- 
nitude below the sensitivity limit of current neutron in- 
strumentation. In this article we continue our program 
of putting this qualitative observation on a quantitative 
footing by converting the neutron cross section to the 
absolute spin susceptibility. 

The optimally doped superconductors YBa 2 Cu307 
and Bi 2 Sr 2 CaCu 2 08+«5 with transition temperatures in 



the 90 K range at first sight appear compatible with a 
weak correlation picture as the normal-state spin suscep- 
tibility is indeed too weak to be measurable with the 
instrumentation currently available |^0|,|ll| , [lq , P2|. How- 
ever, below the superconducting transition temperature 
both compounds develop a remarkably sharp "resonance 
peak" at an energy of 40 meV whose spectral weight is 
comparable to that of antiferromagnetic magnons at the 
same energy p-flOf ■ This discovery, which was followed 
up by related observations in underdoped YBa 2 Cu306+ x 
||12]-|l4| has generated much theoretical work, most of 
which is based on strong-correlation models ]27| |45|] . Ki- 
netic energy effects such as band structure singulari- 
ties 29 , 3(1 El] or the manifestations of interlayer pair 
tunneling along the c-axis |47|| have, however, also been 
implicated. The resonance peak has not been observed 
in conventional superconductors, but some recent data 
on heavy fermion compounds appear to have a similar 
signature pi] . Here we give a comprehensive account of 
this unusual phenomenon that provides important clues 
to the mechanism of high temperature superconductivity. 
We fully describe the evolution of its energy and spectral 
weight with hole concentration and elucidate the relation- 
ship of the resonance peak to the normal state excitations 
in the underdoped regime. 

Our new experiments are generally consistent with pre- 



vious investigations of the spin dynamics in the under- 
doped regime |@-^,0,§1, but due to the combination of 
large single crystals and advanced neutron instrumen- 
tation the energy range could be extended above 100 
meV, comparable to the nearest-neighbor antiferromag- 
netic superexchange energy J in the insulator, and the 
statistics could be substantially improved. Due to the ex- 
tended energy range, we were able to probe not only the 
low energy sector of the spin excitation spectrum which 
evolves out of in- phase precession modes of spins in di- 
rectly adjacent layers of the antiferromagnetic insulator 
YBa2Cu306 (acoustic magnons, Ref. |2|,|52|), but also the 
higher energy sector evolving out of antiphase (optical) 
magnons [ p3[ . In the metallic state we prefer to charac- 
terize these excitations by their symmetry (odd or even, 
respectively) under exchange of the two layers. Observa- 
tion of excitations in both sectors is important because 
it throws light on interlayer interactions within a bilayer 
both in the normal and in the superconducting states, 
an issue that has emerged as an important theme in the 
theory of high temperature superconductivity. 

The article is organized as follows. In Section pi we 
introduce our notation and the absolute scale which un- 
derlies much of the presentation. Technical details of the 
sample characterization and the neutron experiments are 
given in Sections III and [TV|, respectively. Section [V] gives 
a comprehensive account of the results on even and odd 
excitations in the normal state, and the spin gap and 
magnetic resonance peak in the superconducting state. 
The data are discussed in Section VI in the light of re- 
cent model calculations. Preliminary accounts of this 
work have been given in |13,|lq-|lql- 



II. DEFINITIONS AND NOTATION 
A. Magnetic Susceptibilities and Absolute Units 

Throughout this paper, wave vectors are expressed in 
reciprocal lattice units (r.l.u.), that is, Q = Ha* + Kh* + 
Lc* with a* ~ b* ~ 1.63A" 1 and c* - 0.57A -1 . In 
these units, the antiferromagnetic zone center commonly 
referred to as Qo = (7r,7r), is at Q = (|, 5). The two- 
dimensional reduced wave vector, modulo Qo, is denoted 
by q- 

Magnetic neutron scattering probes the imaginary part 
of the magnetic susceptibility tensor |5j| , 
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(i) 



where rg is the classical electron radius, k^ and k/ are 
the wave vectors of the incident and scattered neutron 
with Q = k/ — ki, a, (3 are the spatial components of the 
tensor, /(Q) is the atomic form factor of the orbitals that 



contribute to the inelastic process, IT(Q) is the Debye- 
Waller exponent, and N is the number of spins in the 
system. The susceptibility is related to the spin-spin cor- 
relation function in the form: 

SQ/3(Q,tj) = ^E/* elQ ' (R3 ~ Rs) ~^f^(i)} 






(2) 



where Si is the usual spin operator at site R^. It will 
often be fruitful to refer to this quantity because there is 
a sum rule for it (see Eq. pd| below). 

When there is no magnetic long range order in the 
system and no preferred orientation, the summation in 
Eq. |l| must be rotationally invariant, that is, ^2 a g(3ap — 

^^(q,") = 2Tr( X ^)/3. Therefore, using the 
isotropic susceptibility x" — Tr(%„ /3 )/3, Eq. |l| becomes 
for a magnetically disordered system 
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We had already used this definition previously Jl5|. 
However, in an earlier paper, some of us [[flj used a 
slightly different convention for the dynamical suscep- 
tibility: 



„2 K / 



X 



dVLdE 



= roVr (Q) exp(-2IT(Q))(l + n(u))±-, (4) 



so that x" — (2ttMb/3) x" ■ The current convention 
agrees with the work of Hayden et al. pq | in the metallic 
state of La2- x Sr x Cu04 as well as with their recent work 
in YBa2Cu30g+ x pq ]- Throughout this paper, x"(Q,oj) 
is given in /i^/eV/Junit formula]. 

The signal detected in a neutron scattering experiment 
is proportional to the intrinsic x" convoluted with the 
instrumental resolution function. We extract x"(Qj w ) 
from the data using a model of the resolution function 
and convert the result to absolute units through a calibra- 
tion against selected phonons (see Appendix). Depend- 
ing on the structure of the magnetic spectrum in momen- 
tum space, it is often convenient to further process the 
data by integrating x" over momentum or energy. We 
henceforth often use the spectral weight, here defined as 
J d(Tiui) x"(Qo, w), as well as the Brillouin zone averaged 
susceptibility X2d( w )j a ^ so known as local susceptibility, 
defined as 



X2DM = 



_JAx"(q,«) 



fdPq 



(5) 



The latter quantity is often used for systems with dis- 
persing excitations. In particular, for a two-dimensional 
(2D) Heisenberg antiferromagnet, we have at low energies 
(where the magnon dispersion relation is approximately 
linear in wave vector): 



X^(w)~4S^ x |*!/J~10A«i/eV, 



(6) 



independent of energy. We here used the theoretical 
value of quantum corrections for S — |, Z x = 0.51 [ p7[ , 
although neutron scattering experiments in several un- 
doped cuprates pq, pq| indicate a systematic reduction 
of ~ 30 % of the spin-wave spectral weight as compared 
with the theoretical estimate (possibly due to covalency 
effects). Nevertheless, Eq. ra provides a useful benchmark 
for the "oscillator strength" of normal state magnetic ex- 
citations in the metallic cuprates. The magnetic reso- 
nance peak in the superconducting state, on the other 
hand, is a non-dispersive excitation centered around a 
single point Qo in momentum space. As it is very nar- 
row in energy, the instrumental resolution function car- 
ries out a partial energy-integration so that its spectral 
weight is best expressed in terms of J d(Tuj) x"(Qo, w). 



B. Even and Odd Excitations 

YBa 2 Cu306+ x is a bilayer system, that is, two closely 
spaced copper oxide layers are separated by a much 
larger distance. Information about the magnetic cou- 
pling between the layers has been obtained from studies 
of the spin wave dispersions in insulating YBa 2 Cu3 06.2- 
The superexchange coupling between adjacent bilayers 
is more than four orders of magnitude smaller than the 
primary energy scale, J ~ 100 meV P,|5^]- In metal- 
lic YBa 2 Cu306+ x , where the in-plane correlation length 
is never larger than a few lattice spacings, the bilayers 
are therefore expected to be magnetically decoupled from 
each other, so that the neutron scattering signal is an 
incoherent superposition of the signals arising from indi- 
vidual bilayers. This is in accord with the experimental 
situation. 

The intra-bilayer coupling J± ~ 10 meV |53], on the 
other hand, is only an order of magnitude below the in- 
tralayer superexchange, and it has long been known that 
even in the metallic regime the two layers within a bi- 
layer remain strongly coupled [Q-Q, [JJ • One focus of the 
present study is to further elucidate this magnetic inter- 
layer coupling. 

To this end, we consider an eigenstate \n) of the in- 
plane crystal momentum centered on the layer n. We can 
then define symmetric and antisymmetric combinations 
of states centered on the two layers within a bilayer: 



|*> = (|1> + |2»/V2 
|a) = (|l)-|2))/V2 



(7) 



In this new basis, elementary excitations can be charac- 
terized by transitions between states of the same or op- 
posite symmetry (even or odd excitations, respectively). 
The contribution of these two types of transition between 
spin states to the neutron scattering signal can be varied 
deliberately by adjusting the momentum transfer com- 
ponent Q z perpendicular to the copper oxide layers. In 
particular, with layer states |1) and |2), 



*|l) = d/2|l) 

z\2) = -d/2\2) 
(a, s\e tQzZ \a,s) — cos(Q z d/2) (even) 
(s,a\e lQ ^\a,s) = ism(Q z d/2) (odd) 



(8) 



where d is the intrabilayer distance. Therefore, the inten- 
sities of the odd and even excitations, regardless of their 
in-plane wave vector and energy dependence, follow the 
sm 2 (Q z d/2) and cos 2 (Q z d/2) modulation, respectively, 
and we can write, following Eq. 0, 



X "(Q,o;)=x^ d (q, W )sm 2 (Q,d/2) 

+Xevcn(q^)cOS 2 (Q 2 d/2). 



(9) 



In the antiferromagnetic insulator, Xodd anc ^ Xeven cor ~ 
respond to the acoustic and optical spin-waves, respec- 
tively [|53| . The local and the energy-integrated suscepti- 
bilities as defined in the previous subsection can likewise 
be decomposed into their even and odd components. 



III. CRYSTAL SYNTHESIS AND 
CHARACTERIZATION 

Our samples are three single crystals of volume 
^1 — 3 cm 3 . A precursor powder was made by com- 
busting a stoichiometrically mixed aqueous solution of 
nitrates with sucrose. The powder was then calcined in a 
fully ventilated environment at 940° C for 6 hours to form 
YBCO powders. Crystals were grown using top seeded 
sintering techniques |23[ ] . After fully oxygenating the re- 
sulting crystal by annealing in a pure oxygen environ- 
ment, oxygen was systematically removed by progressive 
annealing under argon flow until the desired fraction of 
oxygen remained. They were then sealed in an evacuated 
quartz tube and annealed for two weeks to ensure homo- 
geneity. We estimate the oxygen content by the weight 
gain/loss of the sample during the anneal, in conjunction 
with comparing to the T c calibration of Cava et al. |39| . 
Several factors contribute to a rather large systematic er- 
ror (Ax ~ ±0.05) in the estimate of the absolute oxygen 
content. For instance, a significant fraction of a foreign 
phase (Y 2 BaCu05, the so-called green-phase) embedded 
in the single crystal matrix makes it difficult to translate 
the weight change of the crystal into a change in stoi- 
chiometry. Three crystals were used: YBa2Cu306.5 with 
T c - 52 K and mass 23.3 g, YBa 2 Cu 3 6 . 7 with T c - 67 
K and mass 14.8 g, and YBa 2 Cu30 6 . 85 with T c ~ 87 K 
and mass 9.5 g. 

The superconducting properties of these samples were 
characterized by SQUID magnetometry. The sharpness 
of the diamagnetic transition is usually a good indicator, 
but by no means a guarantee, of the uniformity of the 
oxygen distribution in a large crystal. A surface sheet of 
high oxygen content material with a high T c may in prin- 
ciple mask more poorly superconducting material in the 
interior. At the expense of reducing the sample volume, 



the SQUID measurements on the YBa2Cu3C>6.5 sample 
were therefore performed on a piece cut from the inte- 
rior of the sample. Details of these measurements have 
been given elsewhere Ju|. On the YBa2Cu30g.7 sample, 
we have, in addition, employed a nondestructive neutron 
depolarization technique that is sensitive to the distribu- 
tion of superconducting transition temperatures through- 
out the bulk of the sample. This method is described in 
Section IV below. 



IV. EXPERIMENTAL DETAILS 

The crystals were aligned such that wave vector trans- 
fers of the form (H , H , L) or (3H, H, L) were accessible 
and loaded into a cryostat. The initial experiments were 
carried out at Brookhaven National Laboratory (BNL), 
on the H4M, H7 and H8 thermal triple axis spectrom- 
eters of the High Flux Beam Reactor. The bulk of the 
work was then carried out on the 2T spectrometer at the 
Laboratoire Leon Brillouin (LLB) for the YBa2Cu306.7 
sample, and on the hot-neutron INI and thermal-neutron 
IN8 spectrometers at the Institut Laue-Langevin (ILL) 
for the YBa2Cu306.5 sample. Additional polarized-beam 
experiments were done on IN20 at the ILL. At BNL, we 
used pyrolytic graphite (PG) (002) monochromators and 
analysers with a PG filter in the scattered beam and neu- 
trons of 30.5 meV fixed final energy. Collimations were 
set at 40-40-80-80. At the ILL, the incident beam was 
monochromated by a vertically focused Cu(lll) crystal 
on IN8, and Cu(200) or Cu(220) with vertical focusing 
were used as monochromators on INI. Vertically and 
horizontally focusing PG(002) was used as analyser. At 
LLB, we used a vertically curved Cu(lll) monochroma- 
tor and a vertically and horizontally focusing PG(002) 
analyser. Heusler(lll) crystals were used in the polar- 
ization work at IN20. No beam collimations were used at 
the ILL and the LLB in order to maximize the benefits 
of focusing. A PG filter was placed behind the sample 
on 2T and IN8, and the final energy was fixed at 14.7, 
30.5, or 35 meV depending on the range of excitation en- 
ergies covered. An Erbium filter was used on INI, and 
the final energy was fixed at 62.6 meV. The energy and 
momentum resolutions varied depending on the instru- 
mental configuration used. Typical values for the en- 
ergy resolution are 1 meV (full width at half maximum, 
FWHM) at energy transfer 5 meV, 5 meV at 50 meV, 
and 12 meV at 100 meV. The vertical resolution was typ- 
ically 0.35 A -1 FWHM. Full four-dimensional resolution 
calculations were used to extract the instrinsic magnetic 
neutron cross section from the raw data. 

The Heusler monochromator used in the polarized 
beam experiments JllLpJ] reflects only neutrons of a spe- 
cific (vertical) spin polarization direction. Before im- 
pinging on the sample, the beam polarization is main- 
tained (vertical field, VF) or rotated by 90° (horizon- 
tal field, HF) by homogeneous guide fields. After scat- 



tering from the sample, the beam polarization is again 
maintained or rotated back by 90°, respectively. The 
beam then traverses a flipper (a set of coils capable of 
flipping the neutron spin polarization by 180°), and the 
final beam polarization is analysed by a Heusler crystal 
which Bragg-reflects only neutrons whose polarization di- 
rection is the same as that of the original beam (after the 
monochromator) . Because of limitations of the apparatus 
the beam polarization is always incomplete and is usually 
parametrized as (FR-1)/(FR+1), where FR is the "flip- 
ping ratio". When the flipper is on, the spin-flip (SF) 
cross section is measured, superposed by a polarization 
"leakage" contribution from non-spin-flip (NSF) scatter- 
ing events (mostly phonon scattering), a contribution 
from nuclear spin incoherent scattering (NSI) , and an ex- 
trinsic background (B). Because of polarization terms in 
the coherent magnetic scattering cross section [M only 
half of the magnetic contribution (M) is measured for 
vertical guide field, whereas for HF the full contribution 
is measured. When the flipper is on and the flipping ratio 
is not too small, one obtains |Q4| 



I 



HF 



M- 



-NSI 



NSF 



Ivf = ~M + ^NSI 



FR 

NSF 



B 



FR 



B 



(10) 



The standard method of extracting the magnetic con- 
tribution to the cross section is to subtract Ivf from Ihf > 
which yields M/2. 

We used this setup on IN20 in order to determine 
the superconducting transition temperature of one of 
our large single crystals (YBa2Cu30g.7) nondestructively 
and in situ. The technique relies on the fact that a 
spin-polarized neutron beam is depolarized when tran- 
siting a region containing a magnetic field whose direc- 
tion varies on a short length scale, such that the neutron 
spins do not follow the change in field direction adiabat- 
ically. Depolarization is negligible only if the parameter 
7] = lB/v^ » 1, where 7 = 1.8 x 10 8 T^sec" 1 is the 
neutron gyromagnetic ratio, B ~ 10 G the applied field, 
v ~ 2000 m/sec the neutron velocity, and 4^ the direc- 
tional variation of B p5[ . The spectrometer is set for a 
nuclear Bragg reflection. The field B is applied at the 
sample position in the same direction as the guide field 
before and after the sample, so that the beam polariza- 
tion is maintained over the neutron flight path and the 
full nuclear Bragg intensity is detected. This remains un- 
changed as the sample is field-cooled through the super- 
conducting transition. Because of the porous microstruc- 
ture of our samples due to inclusions of a second phase, 
Y2BaCu05, magnetic flux penetrates even at low fields 
and gets trapped by microstructural defects such as twin 
boundaries and the Y2BaCuOs inclusions. This flux has 
the same direction as the guide field and therefore does 
not depolarize the beam. The field is subsequently turned 
by 90 degrees at low temperature. Because of flux pin- 
ning, the vortices do not follow this change of field ori- 
entation, so that the net field at the sample position is 



no longer parallel to the guide field. Significant beam 
depolarization results because J§ ~ j^—, where 1 cm is 
a typical sample dimension, and hence ry < 1. The mea- 
sured intensity of the nuclear Bragg reflection is there- 
fore reduced (Fig. |f|). Heating the sample in this state 
and monitoring the Bragg intensity provides a T c -curve 
characteristic of the entire bulk of the sample. The small 
transition width shown in Fig. ffl attests to the high qual- 
ity of the YBa2Cu3C>6.7 sample. 

While neutron scattering with polarization analysis is 
the most powerful tool to study magnetic excitations, 
current instrumentation restricts its use to rather low en- 
ergy transfers (below ~ 50 meV). Moreover, the neutron 
flux is substantially lower than in a standard unpolarizcd- 
beam experiment. For an extensive survey of the mag- 
netic spectrum, and for higher excitation energies, inelas- 
tic neutron scattering experiments without polarization 
analysis are the only option. In unpolarized-beam ex- 
periments, the dynamic structure factor of phonons can 
exhibit various features in momentum space that can be 
mistaken for magnetic fluctuations at the worst, or com- 
plicate background subtraction at the very least. In or- 
der to arrive at a reasonable description of the phonon 
background, we numerically modeled the lattice dynam- 
ics of YBa2Cu30g+x using a simple harmonic interaction 
model. Once the phonon dispersion relations and eigen- 
vectors are known as a function of wave vector, their full 
cross section can be calculated in absolute units. This 
has the added benefit of providing a standard against 
which the magnetic cross section can be calibrated. In 
previous publications, we have demonstrated the success 
of this model in predicting the cross section of a partic- 
ularly important phonon at 42.5 meV (|,|ll[. We have 
since augmented our calibration procedure by including 
several acoustic phonons, and we have further improved 
the parameters used in the simulation. Details are given 
in the Appendix. 

We do not attempt to fit the whole phonon spectrum 
with our somewhat simplified model of the lattice dynam- 
ics which is still not completely understood |3(J. Nev- 
ertheless, the simulation has given us a useful guide to 
extract the magnetic signal in some situations. As an ex- 
ample, Fig. shows a "worst-case" scan at 55 meV, an 
energy range in which the phonon cross section is partic- 
ularly large in the region of momentum space where the 
magnetic cross section is peaked. This excitation energy 
is too large for polarized-beam experiments. The figure 
shows that the model calculation allows us to extract the 
magnetic cross section from the data with some confi- 
dence even in this energy range. The dashed-dotted line 
describes the prediction of the phonon model which ac- 
counts well for the bowl-shaped background of the scan. 
The peak in the center of the scan, on the other hand, is 
not described by the phonon model and can therefore be 
recognized as magnetic. 

Also, in many cases, additional elements such as the 
momentum dependence of the neutron intensity over 



several Brillouin zones, its temperature dependence, its 
behavior upon changing the resolution conditions, etc., 
have been employed to cross-check the determination of 
the magnetic scattering. Such empirical methods have 
been discussed and successfully applied in previous stud- 
ies 0,0. 



V. RESULTS 

Typical constant-energy scans (whose background is 
generally much more benign than that in Fig. 0) are 
shown for YBa2Cu306.5 in Fig. || and for YBa2Cu30g.7 
in Fig. H. Several qualitative features are already ap- 
parent from the raw data. First, while the magnetic 
signal is always peaked at or near Qo = (ir, 71"), the de- 
tailed shape of the profiles in momentum space depends 
on the excitation energy. As the excitation energy in- 
creases, the peak generally broadens, and at least in the 
YBa2Cu306.5 sample it begins to disperse away from Qo 
above 50 meV (Fig. |J). The line shape will be discussed 
in Section |VC| below. Second, the even spin excitations 
(shown in Fig. 0) are fully gapped, as are the even ex- 
citations (optical spin waves) in the antiferromagnetic 
insulator. The data were fitted to single or double Gaus- 
sians, as appropriate, and corrections for the resolution 
function were made. 

The data were further corrected for the magnetic form 
factor of copper. Because of the energy dependence of 
the momentum line shape ( Figs. 0-0), we plot in Figs. Jq 
and U the q-averaged (local) susceptibility (x 2 ' D , Eq. ||) 
of the two samples in the even and odd channels, derived 
from the fitted intensities and widths of the Gaussian pro- 
files. Alternatively, the data can be summarized in terms 
of the peak intensity at Qo, as shown for YBa2Cu30@.7 
in Fig. 0. The local and peak susceptibilities are not 
proportional to each other, because of the energy depen- 
dence of the momentum line shape. Note also that the 
absolute magnitudes of the susceptibilities in Figs. || and 

differ by more than order of magnitude, because the 
susceptibility is always strongly peaked at or near Qo. 
Nevertheless, the qualitative features of Figs. [| and |?] are 
similar. The spin excitations in the YBa2Cu306.85 crys- 
tal are restricted to a much smaller energy range than 
in YBa2Cu306.5 and YBa 2 Cu306.7, as observed previ- 
ously 0. A full spectrum is reported in a forthcoming 
publication p6[. 







A. Magnetic Resonance Peak and Superconducting 
Energy Gap 

One of the most striking features of the data of Figs. 
@a-|7k is a pronounced peak that develops at low tem- 
peratures. Various features of this peak are strongly 
reminiscent of the magnetic resonance peak that was ob- 
served at an excitation energy of 40 meV in optimally 



doped YI^CusQa+x []5| — [Tot] and more recently also in 
Bi2Ca2CuSr208+5 |22|] : it is concentrated around a sin- 
gle point in an energy-wave vector diagram, its intensity 
decreases with increasing temperature, and it occurs in 
the odd channel. This analogy is further bolstered by 
considering the detailed temperature dependence of the 
peak local susceptibility that is shown in Fig. || for both 
samples. Clearly, the intensity is sharply enhanced be- 
low a temperature that, to within the experimental error, 
is identical to T c in both cases. The coupling to super- 
conductivity is another characteristic feature of the reso- 
nance peak in the optimally doped compounds. We can 
therefore clearly associate the enhancement of the dy- 
namical susceptibility in the underdoped samples below 
T c with the magnetic resonance peak. 

As was noted earlier 12 -14], the energy of the reso- 
nance peak decreases systematically as the hole concen- 
tration is lowered from optimal doping. Another differ- 
ence between the optimally doped and underdoped ma- 
terials is that in the latter samples magnetic excitations 
are detectable in the normal state. The buildup of the 
magnetic correlations, though strongly enhanced in the 
superconducting state in a narrower energy window, al- 
ready begins much above T c . Indeed, the response shown 
in Figs. lH^is already broadly peaked in the normal state. 
The norm al state spectra are discussed in detail in Sec- 
tion VB below. 



For now, we focus on the enhancement of the dynami- 
cal susceptibility below T c . In order to elucidate the in- 
fluence of superconductivity on the magnetic excitation 
spectra, we have subtracted the data just above their re- 
spective T c from the data deep in the superconducting 
state. The result is shown in Figs. and HO in the form of 
x"(Qo, w) as well as S(Qo, w). In this context it is fruitful 
to consider the latter quantity because it is constrained 
by a sum rule: 



/ d 3 Q d{hu) Tr(S tt/3 (Q, M) 
JdHJ 



= S{S + 1) per Cu atom 
(11) 



This so-called "total moment sum rule" is strictly valid 
only for the undoped Heisenberg antiferromagnet (with 
spin S = i), and the numerical constant is likely to 
be somewhat reduced as holes are added in the metal- 
lic regime of the phase diagram. Independent of doping, 
however, the total moment sum rule implies that S(Q, w), 
integrated over all momenta and energies, is temperature 
independent. Note, however, that in a metallic system 
the relevant energies extend up to energies comparable 
to the Fermi energy, far beyond the energy range probed 
by our experiment. 

Fig. J9|b shows that for the YBa2Cu306.5 sample the 
enhancement of the spectral weight around 25 meV, 
which we attribute to the magnetic resonance peak, is 
accompanied by a reduction of spectral weight over a 
limited energy range both above and below 25 meV. 
The width of the response in momentum space does not 



change significantly upon entering the superconducting 
state, and the reduction of spectral weight above and 
below the resonance peak compensates the resonant en- 
hancement such that the total moment sum rule is satis- 
fied to within the experimental error. 

The YBa2Cu306.7 sample shows a qualitatively sim- 
ilar behavior, but within the energy range probed by 
our experiment the enhancement of the spectral weight 
around the resonance energy of ~ 33 meV outweighs 
the reduction of spectral weight at other energies. This 
trend continues as the hole concentration in increased. 
In YBa2Cu307, where normal state excitations have not 
been clearly identified, a loss of spectral weight accompa- 
nying the resonance peak has thus far not been observed. 
At least for doping levels exceeding x ~ 0.7, we therefore 
have to postulate a broad and weak continuum of exci- 
tations, perhaps extending up to high energies and not 
directly observed in our experiments, as a "reservoir" of 
quantum states from which the resonance peak is drawn. 

Following these considerations, we have chosen to pa- 
rameterize the spectral weight of the resonance peak by 
the positive component of the difference spectra of Figs. 
I) and [h], f d(huj)Ax+ (Qo, u>), in order to compare dif- 
ferent doping levels. (Note that in the energy range in 
which the enhancement occurs, and over the temperature 
range considered, the Bose population factor 1 + n(cu) in 
Eq. is very close to 1, so that A%" and AS+ are sim- 
ply proportional to each other.) Fig. |ll| summarizes the 
doping dependences of the spectral weight and energy of 
the resonance peak as a function of doping concentra- 
tion. The large error bars attached to A%" in Fig. O 
are in part due to ambiguities arising from the incom- 
mensurate response below the resonance peak which is 
discussed below. While the energy of the resonance peak 
increases with increasing doping level, the absolute spec- 
tral weight is more weakly doping dependent (though a 
clear reduction is observed in the fully oxygenated sam- 
ple YBa2Cu30y). Note also that the dynamical suscep- 
tibility in the normal state decreases strongly with in- 
creasing doping level. A feature not shown in Fig. O 
is the width of the magnetic resonance peak in energy 
which is comparable to the experimental energy resolu- 
tion for YBa2Cu3C>6.7, YBa2Cu3C>6.85 and YBa2Cu3C>7, 
but broadened to an intrinsic width of ~ 10 meV for 
YBa 2 Cu306.5- It was shown in pl[ that a small amount 
of disorder can lead to a drastic broadening of the peak, 
and we cannot exclude disorder in the Cu-0 chain layer 
as the origin of the broadening in YBC06.5- 

Although the q-width of the resonance peak, Aq = 
0.25A , does not depend significantly on doping, it is 
nonetheless interesting to compute the energy and wave 
vector integrated spectral weight of the resonance peak, 
and to compare the result to the total moment sum rule, 
Eq. [ll]. After a two-dimensional Gaussian integration 
of the difference signal, we obtain 0.069/4, 0.056/4, 
0.07/4, and 0.043/4 for YBa 2 Cu 3 6 . 5 , YBa 2 Cu 3 6 .7, 
YBa 2 Cu306.85> an d YBa 2 Cu307, respectively. This 



means that roughly 1-2% of the total magnetic spectral 
weight is redistributed to the magnetic resonance peak 
upon cooling below T c . 

Another important phenomenon at least partly associ- 
ated with the superconducting state is the formation of 
a spin gap in the susceptibility spectrum with decreas- 
ing temperature. In the La2-xSr x CuC>4 superconduc- 
tor, the normal state response is incommensurate and 
a small spin gap (~ 3-6 meV) is observed in the su- 
perconducting state |p3|)|64|. In this system, the size of 
the spin gap seems to be a good indicator of the homo- 
geneity of the samples; impurities apparently introduce 
new states below the spin gap, thereby smearing out the 
superconductivity-induced anomalies in both tempera- 
ture and energy. This has also been confirmed experi- 
mentally in Zn-substituted YBa2Cu30g+ x 65, 36[ . 

We have established the size of the spin gap in the su- 
perconducting state of YBa2Cu3C>6.7 in a polarized-beam 
experiment. Polarization analysis is required for an accu- 
rate determination of the onset of magnetic scattering, as 
optical phonon scattering is rather strong in the vicinity 
of the spin gap. Fig. O shows typical constant-energy 
scans taken with a polarized beam, and Fig. O shows 
the low energy spectrum of YBa2Cu30g.7 derived by fit- 
ting the constant-energy scans to Gaussians and plotting 
the amplitude versus energy. The large spin gap of ~ 17 
meV is testimonial to the quality of the samples. Its size 
agrees well with the measurements of Rossat-Mignod et 
al. on a YBa2Cu30g.69 sample p|. Similar measurements 
for YBa2Cu306.5 yield an energy gap of about 5 meV at 
low temperatures. The doping dependence of the super- 
conducting spin gap is therefore quite different from that 
of the magnetic resonance peak itself, as already found 
previously g-Q]- 



B. Normal State Susceptibility 

Having established the influence of superconductivity 
on the spin excitations, we now describe their develop- 
ment with temperature and doping in the normal state. 
The general trends can be seen in Figs, raando, where the 
susceptibility is given in absolute units up to excitation 
energies of 120 meV. The low-energy spectra in the odd 
channel, which are consistent with previous work |2|-[|,f7|, 
show a broad peak whose amplitude decreases gradu- 
ally with increasing temperature and doping. Figs. pHfl 
demonstrate that the low energy spectral weight is al- 
ready rather small in the normal state, and that the ad- 
ditional depression upon entering the superconducting 
state is rather subtle. (This was established for our sam- 
ples in detailed studies of the temperature dependence 
of the low energies, not shown here.) This is consistent 
with prior neutron observations of the spin pseudo-gap 
in underdoped YBa2Cu30g+x as well as in nuclear 
magnetic resonance (NMR) measurements |6^,g8[. The 
relationship between the magnetic resonance peak and 



this normal state peak will be described in the next sec- 
tion. We here stress that one has to be cautious in inter- 
preting the broad normal-state peak in the odd channel 
simply as a precursor of the magnetic resonance peak in 
the superconducting state. 

Until recently, neutron studies of YBa2Cu306+ x have 
been confined to the odd channel, where excitations can 
be observed at relatively low energies. This situation was 
changed when Reznik et al. [p3| identified even excita- 
tions (optical magnons) in the antiferromagnetic insula- 
tor. Although the coupling between two directly adjacent 
layers, extracted from an analysis of the spin wave dis- 
persions, is only J± ~ 10 meV, antiferromagnetic long 
range order, together with the much larger intralayer ex- 
change coupling of J ~ 100 meV, boosts the gap for 
optical magnons to ~ 67 meV. It is therefore interesting 
to monitor the odd-even splitting in the absence of long 
range order. Figs, pp and pb show that the even spec- 
trum remains fully gapped even in the metallic regime, 
although the gap decreases and broadens with increasing 
doping. 

Similar to the odd channel, the even excitations also 
exhibit a broad peak, albeit centered around a higher en- 
ergy. The temperature evolution of the intensity around 
the peak position is given in Fig. |lj for YBa 2 Cu306.7 
and compared to the intensity in the odd channel at the 
same energy. The even intensity increases markedly with 
decreasing temperature, whereas the odd intensity at this 
energy remains constant to within the error. The marked 
difference in even and odd response functions in this en- 
ergy range is already quite apparent in the data of Fig. |j. 
The gradual increase is reminiscent of the gradual growth 
of the peak intensity in the odd channel in the normal 
state (Fig. |3|), but no anomaly is observed in the even 
channel at T c . 

We summarize the temperature evolution of even and 
odd spectra schematically in Fig. [15|. At high temper- 
atures the dynamical susceptibility is rather featureless. 
Lowering the temperature leads to a gradual enhance- 
ment of the spin susceptibility in both even and odd 
channels at approximately the same rate, but centered 
around different frequencies. Just above T c , both spec- 
tra therefore look very similar, with a broad peak in each 
channel that is merely shifted in frequency. Below T c , the 
parallel evolution of both spectra ceases. The peak in the 
odd channel sharpens abruptly while the intensity in the 
even channel continues its smooth normal-state evolution 
and eventually saturates. 

Finally, we draw attention to a noticeable dip in the 
odd channel excitation spectrum around 55 meV (Figs. |5| 
and o) which we are unable to account for in our phonon 
simulation. Above this minimum the local susceptibil- 
ity rises to a second (less pronounced) maximum. This 
behavior is in qualitative agreement with pulsed neutron 
data in YBa2Cu30g.6 j5q |. (Note, however, that the in- 
tensity above 70 meV reported in Ref. |56| is twice as 
large as ours even though both data sets agree in the low 
energy range). The dip is observed in both underdoped 



samples but is more pronounced in YBa2Cu30g.5. It is 
tempting to associate this feature with the gap in the 
even spectrum that also occurs in this energy range. A 
more detailed analysis of the spectra reveals, however, 
that the gap in the even channel has a stronger doping 
dependence than the 55 meV feature in the odd channel, 
which makes this scenario unlikely. We briefly discuss 
alternative explanations in Section |Vl| below. Unfortu- 
nately, the dip feature occurs in an energy region with a 
nontrivial phonon background (Fig. |2|) and at an energy 
too high for polarization analysis. A detailed experimen- 
tal characterization is therefore difficult. 



C. Momentum Line Shape 

In the context of recent debates about the role of 
"charge stripe" fluctuations in high temperature super- 
conductivity, the momentum line shape of the spin re- 
sponse in YBa2Cu306+ x has emerged as an important is- 
sue. Early indications of an incommensurate response |S|] 
were followed up by more detailed investigations over a 
narrow range of energies and doping levels [|19[ Q ■ Be- 
cause of the coarse vertical resolution of a triple axis 
spectrometer and the nontrivial phonon background, spe- 
cial precautions are required to clearly resolve the in- 
commensurate response in YBa 2 Cu30g +x . This is not 
the primary focus of the present study; ongoing de- 
tailed measurements of this aspect are focused on the 
YBa 2 Cu306.85 sample, where the incommensurate re- 
sponse is most pronounced, and will be reported sepa- 
rately in a forthcoming publication P6| . Here we discuss 
some generic features of the momentum line shape that 
are apparent without high momentum resolution. 

Independent of the detailed line shape, an important 
parameter for the purposes of the present study is the 
overall extent in momentum space of the magnetic re- 
sponse at a given energy. For instance, this quantity en- 
ters into our determination of the local susceptibility of 
Figs. @ and O. This overall width was extracted from the 
constant-energy profiles by fitting them to either a single 
Gaussian centered at Qo = (jr, 7r), or if necessary, to two 
Gaussians symmetrically displaced around Qo- Fig. |lq 
gives the result of this procedure for our YBa2Cu306.7 
sample. While most profiles are adequately fit by a broad 
single peak, we paid particular attention to the energy 
range around 24 meV, just below the resonance peak , 
and were able to confirm the recent observation [19] of 
an incommensurate response. The weakly energy depen- 
dent q- width below the resonance peak has been noticed 
before [0,fln4,p2| and may reflect an incommensurate re- 
sponse which is unresolved due to the coarse resolution. 
At the resonance energy, the overall q-width goes through 
a minimum which is particularly noticeable when consid- 
ering the resonant response only (that is, the q-width of 
the additional intensity peak below T c , Fig. |10|) p0| . 

The pronounced intrinsic broadening at high energies 



above 50 meV cannot be attributed to resolution ef- 
fects and is reminiscent of the dispersion-like behavior 
observed in YBa2Cu306.5 (Fig. ||, see also p5|) which 
in turn resembles the dispersion of antiferromagnetic 
magnons. It is interesting to notice that the onset of this 
broadening coincides with the peak-dip feature in Figs. 
ra and o. In this high energy region, the primary effect of 
hole doping therefore appears to be a strong broadening 
of the magnetic response. Overall, Fig. H9 is in good 
agreement with the recent report by Arai et al. |6l| on 
a sample with a similar T c . (Note, however, that these 
authors chose to fit all of their constant-energy profiles 
to double peaks.) 



VI. DISCUSSION 

The most striking feature in the neutron spectra of 
YBa 2 Cu306+ x is the magnetic resonance peak whose ex- 
perimental properties are summarized in Fig. O. This 
phenomenon has been addressed in numerous theoreti- 
cal studies over the past few years. One purpose of the 
present article is to present a summary of the current 
state of experimental information on the resonance peak. 
The other purpose is to stimulate further theoretical work 
on the interplay between normal state excitations and the 
resonance peak in the underdoped regime. We discuss 
these two aspects in turn. 

Compared to the broad continuum expected for ordi- 
nary metals, the resonance peak involves an extraordi- 
narily small volume of phase space, d-wave supercon- 
ductivity helps explain this observation: The coherence 
factor in the neutron scattering cross section has a pro- 
nounced maximum at a wave vector connecting two lobes 
of the d-wave gap function with opposite sign of the or- 
der parameter H. (An early interpretation based on a 
different gap symmetry p7| has been ruled out by a va- 
riety of other measurements [p9[.) It has since turned 
out that d-wave superconductivity alone is not sufficient 
to explain the sharpness and spectral weight of the res- 
onance peak, and that other factors further narrow the 
phase space involved in this process. Strong Coulomb 
interactions between electrons are the most likely factor. 

Some of the features of the resonance peak can be 
described in models that do not include such interac- 
tions. For instance, early photoemission studies have 
suggested an extended saddle point singularity near the 
Fermi level J70| . It has been shown that this band struc- 
ture anomaly leads to a peak in the joint density of states, 
which in turn produces a peak in the unrenormalized sus- 
ceptibility x'o pOJpill- The peak grows as the gap opens 
in the superconducting state. The peculiar gap struc- 
ture in the interlayer tunneling model of high tempera- 
ture superconductivity also leads to a rather sharp peak 
in x'o m the superconducting state p7| . The peak posi- 
tion is determined by a combination of the energy gap 
and the chemical potential (see, e.g., Ref. p3J) so that 



even in this minimalist model the resonance peak and 
the superconducting energygap do not have to show the 
same doping dependence |56|. The resonance peak en- 
ergy decreases as the hole concentration is lowered from 
the optimally doped state (Fig. O) while the energy gap, 
which is directly determined by photocmission and other 
techniques, remains constant or increases (see, e.g., JnJ). 
However, the interplay between the pseudo-gap measured 
in the normal state and the superconducting gap is not 
well understood at present in these experiments. 

There are, however, strong indications that electron 
correlations are necessary for an adequate description of 
the resonance peak. First, the absolute magnitude of 
the dynamical susceptibility, determined experimentally 
in Rcfs. (ll],[15| , is too large to be consistent with the rela- 
tively subtle enhancement assumed in Refs. 1 46,47]]. This 
does not necessarily rule out the mechanisms proposed in 
this work, but it strongly suggests that Coulomb corre- 
lations are also needed to reproduce the experimentally 
measured x" ■ Second, the sharpness of the resonance 
peak close to optimal doping indicates a true collective 
mode. This is indeed the direction that most of the recent 
theoretical work has taken. 

The theoretical description of this collective mode re- 
mains a subject of considerable controversy. In one the- 
oretical approach, the peak is identified with a resonance 
in the particle-particle channel. Proponents of this model 
JhJ, later embedded into a more comprehensive theory 
of high temperature superconductivity EILE2|, early on 
predicted the resonance spectral weight to within a fac- 
tor of two of the measured number for optimally doped 
YBa2Cu306+ x ) PI - In the framework of this theory, the 
resonance can be viewed as a pseudo-Goldstone boson 
of a new symmetry group encompassing antiferromag- 
netism and d- wave superconductivity. The symmetry 
becomes exact at the quantum critical point separating 
both phases, and the energy of the Goldstone boson ap- 
proaches zero in accord with the data. This scenario is 
also consistent with the tendency of the resonance spec- 
tral weight to increase with decreasing hole concentra- 
tion (Fig. [fl]). All of these features are, however, also 
consistent with a more conventional description in which 
the resonance occurs in the particle-hole channel. In this 
approach, the resonance is identified with a spin exci- 
ton that is stabilized when decay channels are removed 
in the superconducting state. Recent theoretical studies 
have shown that the evolution of the resonance energy 
and spectral weight with doping (Fig. [ll]) as well as the 
broadening of the resonance in the deeply undcrdoped 
regime can be explained in this framework 36 , 4q|. It 
has also emerged that a remarkably consistent picture of 
ARPES and neutron data can be obtained in this frame- 
work p6|,|78| , [7^ ]. An interesting suggestion put forth in 
this context is to interpret the weakly doping dependent 
dip-peak feature of Figs, ra and pks a signature of the 
bare superconducting energy gap, while the resonance 
peak is pulled below the gap by interactions [f79| . Further 
work is necessary to establish whether the onset of the 



mo mentu m-space broadening in this energy range (Sec- 
tion VC ) can also be explained in this model. Finally, 
it has been proposed that the resonance peak originates 
from localized electrons in domain walls between charge 
stripes, separate from those that form the superconduct- 
ing condensate p8[ . While the large spectral weight of 
the peak finds a natural explanation in this approach, it 
is much harder to account for its sharpness in both energy 
and momentum and its strong coupling to superconduc- 
tivity. Further, the evidence of charge stripes is weak in 
optimally doped YBa 2 Cu307 where the resonance peak 
is most pronounced. 

We do not wish to reiterate the purely theoretical ar- 
guments that were advanced in this debate (see, e.g., 
lH, H,[|| H,H)- Rather, we point to the additional 
experimental information presented in this article, much 
of which still awaits a theoretical explanation. We first 
disc uss t he normal state susceptibility described in Scc- 
VB and summarized schematically in Fig. Hi 
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broad peak in the odd channel that grows with decreasing 
temperature (already been observed in early work on un- 
derdoped YBa 2 Cu306+ x d~§|;0) signifies a new energy 
scale in the normal-state spin excitation spectrum, dif- 
ferent from the superexchange interaction J that sets the 
energy scale in the insulator. It was pointed out ]37,75 
that damping of antiferromagnetic spin waves (seen as 
long-lived excitations in the insulator) introduces such a 
new energy scale, namely ~ Ja/t; (where £ ~ 2a — 5a, 
depending on doping, is the spin-spin correlation length), 
above which the excitations are not substantially differ- 
ent from those of the insulator. This approach can de- 
scribe some aspects of the data, especially the dispersion- 
like behavior observed at high energies in YBa2Cu3C>6.5. 
The response in the even channel is also qualitatively con- 
sistent with a description based on damped spin waves. 
When £ is short, one may expect the even and odd re- 
sponse functions to be shifted by J±£/a, the energy cost 
for flipping a correlated patch of spins in a single layer 
while keeping the spins in the neighboring layer fixed. 
(This simple estimate was confirmed by numerical sim- 
ulations on finite-sized systems.) The actual situation 
(Figs. H and |6|) interpolates between this limit and the 
antiferromagnetic long range ordered state with an opti- 
cal magnon gap of 2*Jj~J±_ ~ 67 meV j|j|. 

The temperature dependence of the broad normal- 
state peak in the pseudo-gap regime and its relation to 
the magnetic resonance peak are more difficult to de- 
scribe. Of course, short range dynamic spin correlations 
centered around Qo = (tt, tt) are expected in a variety 
of microscopic strong-correlation models. The challenge 
is now to use these models to make detailed predictions 
of the temperature evolution of the magnetic spectra in 
both the superconducting and the normal states. Phe- 
nomenologically, the low-energy "overdamped spin wave" 
response grows with decreasing temperature. The broad 
peaks in both even and odd channels grow gradually and 
approximately at the same rate, while at the same time 
a spin pseudo-gap opens in the odd channel as already 



revealed by earlier neutron [g, [7| and NMR |67]] experi- 
ments. Although it is hard to associate a characteristic 
temperature with this crossover phenomenon, our data 
are certainly consistent with models that predict a cor- 
related spin liquid phase below a doping dependent "co- 
herence temperature" , often termed T* |7j| . 

It is important to draw a distinction between the grad- 
ual growth of the broad normal state peaks in both 
even and odd channels on the one hand, and the abrupt 
increase of the intensity in the superconducting state 
that occurs in the odd channel only on the other hand. 
This is seen most clearly in the high-statistics data of 
Fig. ||a. Though both are clearly related, a recent re- 
interpretation of the previously observed broad normal- 
state peak in the odd channel as a "pre-transitional pre- 
cursor" of the resonance peak is therefore too simple- 
minded [p6[ . There is, moreover, no basis for a sepa- 
ration of resonant and non-resonant response functions 
in the normal state as proposed in Ref. J5(J. Promis- 
ing steps in the direction of a microscopic description of 
the spin dynamics in both normal and superconducting 
states have recently been taken (see, e.g., Ref. p3J)- 

This also bears directly on the issue of the origin of 
the superconducting condensation energy that has re- 
cently moved into the foreground of the high-T c debate. 
Scalapino and White |72| derived a formula that relates 
the exchange energy (associated with the J-term in the 
t — J model) to the magnetic excitation spectrum. If 
the complete spectrum is known in absolute units, the 
change in exchange energy between normal and super- 
conducting states can be evaluated and quantitatively 
compared to the superconducting condensation energy 
determined in specific heat experiments |77y . In opti- 
mally doped YBa2Cu3 06+ x , the resonance peak is the 
only experimentally discernible feature in the magnetic 
spectrum, and its spectral weight is known in absolute 
units |llj| . Since the normal state spectrum is unknown 
except for an upper bound guaranteeing that the ampli- 
tude of the normal state excitations is significantly be- 
low the resonance amplitude [|I0[ pTl|, jig] , it is reasonable 
to apply the Scalapino- White formula to the resonance 
peak only. This is the approach adopted by Demler and 
Zhang |£3[ who found that the magnetic energy stored in 
the resonance peak equals the superconducting conden- 
sation energy to within an order of magnitude. Temper- 
ature dependent changes in other parts of the spectrum, 
though not visible in the neutron experiments, could of 
course lead to quantitative modifications in this analy- 
sis. The full implications of this analysis for the pairing 
mechanism of high temperature superconductivity are, 
however, still under debate. Others J74J have argued, 
for example, that the neutron peak provides a measure 
of the condensate fraction, rather than the condensation 
energy, of the superconducting state. 

The experimental situation is much less straightfor- 
ward in the underdoped regime. We have shown here 
that the amplitude of the magnetic response in the nor- 
mal state is at least as large as the resonance amplitude, 



and that it is spread over a much larger energy range. If 
one wants to extend the Demler/Zhang analysis into the 
normal state of underdoped YBa2Cu306+ x , it is there- 
fore erroneous to focus exclusively on the energy at which 
the resonance eventually develops in the superconduct- 
ing state | p6| , especially since a separation of the normal- 
state response into resonant and nonresonant contribu- 
tions is entirely arbitrary. It is also inappropriate to ne- 
glect the optical channel in such an analysis, because its 
amplitude and temperature dependence in the normal 
state are at least equal to those in the superconduct- 
ing state. Because the Scalapino- White formula differs 
from the total moment sum rule (Eq. O) only through a 
momentum-dependent form factor, any gain in magnetic 
energy must result from temperature dependent changes 
of the momentum line shape which are not even consid- 
ered in Ref. |pq| . It is easy to see that the outcome of 
an adequate analysis of the magnetic energy stored in 
the normal state spin excitations and its relation to the 
specific heat would likely be qualitatively different from 
the one given in Ref. |56| . In particular, the normal- 
state electronic specific heat in the pseudogap regime of 
the YBa2Cri306+ x phase diagram is actually lower than 
at optimum doping fTn, whereas the magnetic spectral 
weight and its temperature dependence are clearly larger 
in underdoped samples. Because of insufficiently accu- 
rate information on the temperature dependence of the 
momentum line shapes, we have not used our spectra to 
attempt such an analysis ourselves. 

In summary, we have presented a comprehensive ac- 
count of the magnetic resonance peak and the normal 
state spin excitations in underdoped YBa2Cu30g+ x . We 
hope that the detailed description of the magnetic spec- 
tra presented here will provide an improved basis for 
models of electronic correlations in the cuprates. 
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APPENDIX A: PHONON SIMULATION 

Empirical descriptions of the eigenvibrations in a sin- 
gle crystal are instrumental in obtaining the absolute 
magnitude of the dynamical susceptibility in our mea- 
surements. However, as there are 39 phonon branches 
in YBa2Cu307 spread over an energy range of roughly 
70 meV, individual phonon branches are usually diffi- 
cult to resolve (with a typical energy resolution of 5 
meV). Fortunately, in past successful experiments, we 
have isolated one particularly strong phonon branch at 
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42.5 meV |9|,|llJ. Since then we have further reinforced 
confidence in our model by identifying and simulating 
acoustic phonon branches. Acoustic phonons are attrac- 
tive for absolute unit measurements because their inten- 
sities are strongly related to the Bragg peak intensities. 
Unlike optical phonons, the most important quantity for 
the acoustic phonons is just one number-the limiting 
value of the structure factor as it approaches the Bragg 
condition. The convolution of their "cone shape" dis- 
persion curve with the instrumental resolution has been 
performed analytically as well as by using an efficient 
sampling technique based on the fast Fourier transform 
to perform a four dimensional numerical convolution [ p0[ . 
We have extensively used a longitudinal acoustic phonon 
branch around the (006) Bragg condition. Our simula- 
tion gives a structure factor of 4.3 x 10 _18 c 2 eV _1 which 
is consistent with our earlier optical phonon calibration 
measurements. We further point out that the calibration 
using phonon scattering from the sample avoids the dif- 
ficulties related to existence of "green phase" impurities 
as well as inclusions of smaller crystallites inevitable in 
these large YBa2Cu30g+ x samples. This may contribute 
to systematic errors when calibrating against a vanadium 
standard |55|.Kg|. 
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FIG. 1. T c curve of the YBa2Cu306.7 sample using the 
depolarization technique as described in the text. The spec- 
trometer is set for the (006) nuclear Bragg reflection. The 
flipping ratio, defined as the ratio of neutrons scattered with- 
out spin flip to those that are scattered with a spin flip, is in- 
finite at a nuclear Bragg reflection in an ideal setup, but finite 
(here: 16) because of incomplete beam polarization. With the 
sample field cooled into the superconducting state and then 
rotated by 90°, the nipping ratio is reduced by trapped flux. 



FIG. 2. A scan of even channel excitations at 55 meV 
that exemplifies the usefulness of numerical simulations of the 
lattice dynamics. The dotted line, derived from the phonon 
simulation program described in the text, provides a good 
description of the background features. The peak in the center 
is magnetic. The solid line is a guide-to-the-eye. 



FIG. 3. Typical constant-energy scans taken on the 

YBa2Cu30e.5 sample in the odd channel. The solid lines are 
the results of fits to Gaussian profiles. The energies of the data 
in the three panels are evenly spaced (50 meV, 65 meV, and 80 
meV). The dotted lines that join the split peaks go through 
(H — 0.5r.l.u., E = OmeV) and illustrate a spin- wave- like 
dispersion, as described in the text. 
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FIG. 4. Typical constant-energy scans obtained on the 
YBa2Cu3C>6.7 sample, in (a) odd and (b) even excitation 
channels. The top panel of (b) shows ILL data, the rest are 
from LLB experiments. The lines are the results of fits to 
Gaussian profiles. 



FIG. 5. Local (2D wave vector averaged) susceptibility of 
YBa2Cu306.5 in the (a) odd and (b) even excitation channels, 
obtained by fitting constant-energy scans.. Solid lines are 
guides-to-the-eye. The dotted line indicates the low- energy 
local susceptibility of the 2D Heisenberg antiferromagnet (Eq. 



|) 



FIG. 6. Local (2D wave vector averaged) susceptibility of 
YBa2Cu306.7 in the (a) odd and (b) even excitation channels. 
Solid lines are guides- to-the-eye. The dashed line superposed 
on the top (12 K) spectrum reproduces the solid line for the 
70K data and illustrates the effect of entering the supercon- 
ducting state. The dotted line indicates the low-energy local 
susceptibility of the 2D Heisenberg antiferromagnet (Eq. ph. 
(from j§). 

FIG. 7. Same data as in Fig. o, but showing the intensity 
at Qo = (it, 7r) instead of the local susceptibility. 



FIG. 12. Typical constant-energy scan at 24 meV 

taken with polarized neutrons in the spin-flip channel on 
YBa2Cu30e.7 in (a) the superconducting and (b) the normal 
state. Open circles are horizontal field (HF) measurements 
taken after the sample was field cooled with the field paral- 
lel to the scattering vector. Solid circles are measurements 
taken after the sample was cooled in vertical field (VF), per- 
pendicular to the scattering plane. For magnetic scattering, 
the background-corrected HF intensity should be twice the 
VF intensity because of polarization factors in the magnetic 
neutron cross section (Eq. hoh . The solid lines indicate that 
this is indeed the case. 



FIG. 13. Determination of the spin gap in YBa2Cu306.7 
by polarized neutrons. The lines are guides-to-the-eye. Open 
square symbols are from direct subtraction between HF and 
VF intensity as in the previous figure. Solid square symbols 
are HF polarized beam background from curve fitting such 
data. Open circles are the amplitudes resulting from these 
fits. Note that the scales on left and right vertical axis are 
the same. 



FIG. 14. Local susceptibility at 65 meV in YBa2Cu306.7 
in both channels, showing their different temperature depen- 
dences. Note the absence of any anomaly at T c . Lines are 
guides-to-the-eye. 



FIG. 8. Temperature dependence of the local susceptibil- 
ity at the energies at which the magnetic resonance occurs in 
both underdoped samples. Superconducting critical tempera- 
tures are marked with arrows. The solid symbols are obtained 
using polarized neutron scattering techniques. 



FIG. 9. Difference of (a) the susceptibility x"(Q0)W) and 
(b) the spin-spin correlation function, S(Qo,w) in the super- 
conducting and normal states, for YBa2Cu306.5. Normal 
state data are taken at 60 K, just above T c = 52 K. 



FIG. 10. Difference of (a) the susceptibility x"(Qo,<^) 
and (b) the spin-spin correlation function, S(Qo,u>) in the su- 
perconducting and normal states, for YBa2Cu30e.7. Normal 
state data are taken at 80 K, just above T c = 67 K. The open 
circles represent data taken at BNL, the filled circles are data 
taken at LLB. 



FIG. 15. Schematic diagram summarizing the tempera- 
ture evolution of the magnetic excitation spectra of the un- 
derdoped YBa2Cu306+ x compounds (a) in the normal state, 
and (b) across the superconducting transition. The even ex- 
citations evolve smoothly with temperature, parallel to (a). 
The odd channel excitations undergo an abrupt sharpening 
at the resonance energy across T c . 



FIG. 16. Overall momentum width of the magnetic re- 
sponse i n YB a2Cu3Q6.7, obtained as described in the text 
(Section VC). Particular attention was paid to the con- 
stant-energy profile at 24 meV, where an incommensurate 
response was identified as first observed by Dai et al. ]19| , 
Note the pronounced broadening above 50 meV. 



FIG. 11. A synopsis of (a) the superconducting transition 
temperature, (b) the energy-integrated spectral weight eval- 
uated at Qo = (77,77), and (c) of the energy of the magnetic 
resonance peak in the two underdoped samples, compared to 
those of the optimally doped sample with T c = 93 K. Hori- 
zontal error bars indicate a conservative confidence level for 
the oxygen content. 
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